home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dlanv2.z / dlanv2
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAANNNNVVVV2222((((3333FFFF))))                                                          DDDDLLLLAAAANNNNVVVV2222((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLANV2 - compute the Schur factorization of a real 2-by-2 nonsymmetric
  10.      matrix in standard form
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )
  14.  
  15.          DOUBLE         PRECISION A, B, C, CS, D, RT1I, RT1R, RT2I, RT2R, SN
  16.  
  17. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  18.      DLANV2 computes the Schur factorization of a real 2-by-2 nonsymmetric
  19.      matrix in standard form:
  20.  
  21.           [ A  B ] = [ CS -SN ] [ AA  BB ] [ CS  SN ]
  22.           [ C  D ]   [ SN  CS ] [ CC  DD ] [-SN  CS ]
  23.  
  24.      where either
  25.      1) CC = 0 so that AA and DD are real eigenvalues of the matrix, or 2) AA
  26.      = DD and BB*CC < 0, so that AA + or - sqrt(BB*CC) are complex conjugate
  27.      eigenvalues.
  28.  
  29.  
  30. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  31.      A       (input/output) DOUBLE PRECISION
  32.              B       (input/output) DOUBLE PRECISION C       (input/output)
  33.              DOUBLE PRECISION D       (input/output) DOUBLE PRECISION On
  34.              entry, the elements of the input matrix.  On exit, they are
  35.              overwritten by the elements of the standardised Schur form.
  36.  
  37.      RT1R    (output) DOUBLE PRECISION
  38.              RT1I    (output) DOUBLE PRECISION RT2R    (output) DOUBLE
  39.              PRECISION RT2I    (output) DOUBLE PRECISION The real and
  40.              imaginary parts of the eigenvalues. If the eigenvalues are both
  41.              real, abs(RT1R) >= abs(RT2R); if the eigenvalues are a complex
  42.              conjugate pair, RT1I > 0.
  43.  
  44.      CS      (output) DOUBLE PRECISION
  45.              SN      (output) DOUBLE PRECISION Parameters of the rotation
  46.              matrix.
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.